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Abstract 

The flux as measured by the Josephson effect in a SQUID-like con- 
figuration with a ferromagnetic core inserted into its center, is shown to 
be sensitive to the vector potential arising from the central ferromagnetic 
core, even when the core is covered with a superconducting material that 
prevents any magnetic field lines from ever reaching the perimeter of the 
SQUID-like configuration. This leads to a macroscopic, Aharonov-Bohm- 
like effect that is observable in an asymmetric hysteresis loop in the re- 
sponse of the SQUID-like configuration to an externally applied magnetic 
field. 

In our ongoing experiments pQ, we are using a ferromagnetic torus 
coated with a superconductor (lead) £2] , which is topologically linked with 
another superconducing (SC) ring made out of a different material (alu- 
minum), which has a lower critical temperature. However, to understand 
this experiment more easily, let us first consider the simpler cylindrically- 
symmetric configuration of Figure I, which depicts a straight, infinitely 
long solenoidal configuration consisting of a cylindrical ferromagnetic core 
(yellow) possessing a permanent magnetization M, which is covered with 
a thick coating of SC lead (blue), in order to prevent, due to the Meissner 
effect, any stray magnetic flux lines from escaping from the core. This 
core in turn is placed at the center of a thin cylindrical shell of SC alu- 
minum (red), which has a lower critical temperature than that of SC lead 
(blue). 

Let us further assume that the temperature has been cooled to near 
absolute zero, and that we are slowly applying a uniform magnetic field 
H onto the entire assembly starting from zero field. The magnetostatic 
flux filling the ferromagnetic core due to its permanent magnetization will 
give rise to a static contribution to the circular lines of the vector poten- 
tial A, as illustrated by the dashed circle in Figure 1. However, note 
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Figure 1: A solenoidal configuration consists of an infinitely long ferromagnetic 
core (yellow) coated by a thick layer of lead (blue), at the center of a thin 
cylindrical shell of aluminum (red). H is an externally applied magnetic field, 
M the magnetization of the ferromagnetic core, and A the vector potential. 
As the temperature is steadily lowered, the lead (blue) coating first becomes 
superconducting, and totally confines the ferromagnetic flux to the interior of the 
core via the Meissner effect, so that the aluminum ring detects no B field arising 
from the core, but detects the A field. As the temperature is further lowered, 
the aluminum (red) ring becomes superconducting, and acquires a persistent 
current which measures the core's A field via a SQUID-like flux trapping effect. 
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that due to the Meissner effect arising from the thick lead coating, the 
magnetostatic field B arising from the ferromagnetic core vanishes in the 
spatial region in which the cylindrical shell of aluminum material resides, 
so that the electrons in this shell never experience any classical forces. 
This is true both above and below the critical temperature of the alu- 
minum shell. Nevertheless, like in the Aharonov-Bohm effect, there exists 
a nonlocal quantum influence of the electrons within the ferromagnetic 
core upon the electrons within the superconducting shell, since persistent 
currents within the core can induce persistent currents within the shell 
that are proportional to the size of the vector potential arising from the 
ferromagnetic core. 

The probability current density in quantum mechanics for a single 
neutral particle of mass m moving through the vacuum is given by 



= ^ (V^Pop^) + cc. (1) 

where %p = %p (r, i) is the wavefunction evaluated at the spacetime point 
(r, t), and where the momentum operator in the configuration space rep- 
resentation is 

Pop = -V (2) 

i 

From the time-dependent Schrodinger equation, the continuity equation 

v.j + | (S, 

follows, in which 

p = V*V = V*(r,t)V'(r»t) (4) 
is the probability density for finding the particle located at (r, t) upon 
measurement. The continuity equation <[3j means that probability is con- 
served. 

Let us generalize the above expression fTJ for the probability current 
density j (r, t) to that of a particle with charge q moving through the 
vacuum in the presence of a vector potential A (r, t), by applying the 
minimal coupling rule [T] 

Pop — > Pop - qA (5) 
to JTJ). Then one finds that for the particle with charge q 

This expression for the probability current density j (r, t) for the charge 
q moving through the vacuum generalizes to an expression for the super- 
current density J s (r, i) for Cooper pairs moving through a SC, which can 
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be obtained from the complex order parameter ^ (r, t) of the Ginzburg- 
Landau theory as follows: 

J ' = 2^H^-* A )*) +C - C - W 

where q — 2e is the charge and m = 2m e is the mass of the Cooper pair. 

Next, let us express the complex order parameter in polar form as 
follows: 

* = VP^ (8) 

where 

p = *** (9) 

is the probability density of Cooper pairs, and where cf> is the local phase 
of the Cooper pair's "macroscopic wavefunction" \&. 

Since the local positive charge density of the background ionic lattice 
must be exactly locally compensated by the local negative charge density 
of the Cooper pairs in order to preserve the overall charge neutrality of 
the SC material (which must be true under all conceivable experimental 
circumstances), we require that 

\p\ = IPioniciatticd = constant (10) 

independent of space and time. Therefore, to an extremely good approx- 
imation, the density of Cooper pairs 

p{T,t)=po (11) 

is independent of (r,t); only their phase <f) = cj>(r,i) is allowed to vary 
with space and time inside the SC. Hence the polar decomposition ((8} 
becomes 

*{T,t) = y /p5e i +™ (12) 
Therefore the supercurrent density ([7| becomes 

Js = i (^WT^ (~ (iV0) - gA) y/pSe* 4 } +c.c. 

Q (iV0) - qAj 



go ffi 
2m 



+ c.c. 



= £1 ( ftV - qA) = p v s (13) 
m 

where v 3 is the superfluid velocity of the Cooper pairs. Hence the super- 
current density becomes 



Po 

Js -COVs — 

and therefore 



J s = p oVs = ^ {HV<t> - qA) (14) 
m 



rnvs = hV4> — qA (15) 

where rav, is the kinetic momentum, h\7(f> is the canonical momentum, 
and qA is the electromagnetic momentum of the Cooper pairs. 

The concept of fluxoid quantization follows from (|15[l , by integration 
over a circular path inside the material of the red cylindrical shell of Figure 
1 as follows: 

, ■ dl = h & X7<j) ■ dl - q & A ■ dl (16) 
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Due to the single- valuedness ol complex order parameter ^, it follows that 



• dl = 2im, (n = 0, ±1, ±2, ...) (17) 
and therefore that 

q <j) A ■ dl+m <j) v s ■ dl = 2imh = nh, (n = 0, ±1, ±2, ...) (18) 
Let us define the fluxoid as follows: 

<l"fluxoid = f (A+— v s ) dl = (19) 



$ = ib A • dl (20) 



where 



is the flux, and where 

s ■ dl (21) 

is the circulation of the superfluid, which can become important when 
the cylindrical shell becomes thin compared to the penetration depth of 
the SC (here, aluminum), as happens when the applied magnetic H field 
approaches the critical field H c . It follows from (|18[) that the fluxoid obeys 
the following quantization condition: 

$fluxoid = n— =n$ ,(n = 0,±l,±2,...) (22) 

where 

$o = y = 2.07 x 10~ 15 webers (23) 

is the quantum of flux. However, if the cylindrical shell is thick compared 
to the penetration depth, then the circulation k for the superfluid deep 
inside the red ring of Figure 1 becomes negligible, and fluxoid quantization 
reduces to the usual flux quantization condition, where 

$ flux = I A ■ dl = $„ = = n$ , (n = 0, ±1, ±2, ...) (24) 

Figure 2(a) depicts a SQUID-like model for the nonlinearity that oc- 
curs within the red cylindrical shell in Figure 1 as the externally applied 
magnetic field passes through the effective critical field of the SC mate- 
rial (i.e., the aluminum ring plus an effective Josephson junction). This 
nonlinearity is essential for a breakdown of the superposition principle for 
electromagnetic fields in a vacuum, which permits the detection of the 
vector potential arising from the ferromagnetic core of the configuration 
shown in Figure 1. Note that we have replaced the red cylindrical shell of 
Figure 1 by a red ring in the form of a thin, flat circular gasket which has 
a gap in it, as shown in Figure 2(b), so as to be able to incorporate the 
Meissner-shielded ferromagnetic core into its center during the assembly 
process. 
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Figure 2: (a) SQUID-likc model for the production of a persistent supercurrent I 
within a SC ring (in red) by a total enclosed flux $, when the ring is interrupted 
by a Josephson junction possessing a critical current Jj. The self-inductance L of 
the ring will lead to the onset of a remnant field and hysteresis, (b) Expanded 
view of the Josephson junction. An aluminum plate (pink) with an oxidized 
sapphire surface (AI2O3), when screwed tightly onto the gapped aluminum ring 
(i.e., in the form of a flat gasket with a gap in it), makes the ring into a closed 
SC circuit containing a single, effective Josephson junction in it. 
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One justification for using this SQUID-like model [3] [4] [5] goes as fol- 
lows: The ring gasket of aluminum in Figure 2(a) must have a large 
enough gap cut into it in order to enable the ferromagnetic torus to slip 
by the gap and enter into the center of the SC ring, in order to realize 
the configuration of Figure 1. But in order to form a complete SC cir- 
cuit, the gap in the aluminum gasket will be filled in by pressing onto the 
top of the gap an overlapping rectangular aluminum plate (shown in pink 
in Figure 2(b)), which covers the gap. Thus a topological "link" can be 
formed between a "ring" and a "torus" pQ, with the same topology as that 
in Figure 1. 

Two small screws, which are drilled into the two overlap areas between 
the rectangular plate and the gapped gasket, will be used to press tightly 
the flat, rectangular aluminum plate onto the gapped region of the flat 
aluminum gasket, and thus to complete a SC circuit. However, since the 
aluminum material of the ring and of the plate will naturally form an 
oxide (AI2O3) layer on its surfaces, two oxide-layer Josephson junctions 
will naturally be formed at the two overlap areas between the plate and 
the ring, after the rectangular plate has been pressed tightly onto the 
gapped ring by means of these two small screws. The resulting two in- 
series, oxide-barrier Josephson junctions of Figure 2(b) can be thought of 
as being a single, effective Josephson junction interrupting the ring. This 
effective Josephson junction will obey the following sinusoidal relationship 
®: 

I = 7jsin0 (25) 

where / is the supercurrent flowing through the junction, Ij is the effec- 
tive Josephson critical current, and <j> is the phase difference across the 
junction, which, for configuration shown in Figure 2(a), is given by the 
Aharonov-Bohm phase 

(j> = ~ d> A • dl = 2tt-5- oc $ (26) 
h J $0 

where $ is the total flux enclosed by the ring and 

*o = A (27) 

is the quantum of flux. Thus the phase shift <f> is directly proportional to 
the total flux <3? enclosed by the ring. Note that this implies that there 
exists a nonlinear relationship exists between the supercurrent / and the 
total flux on account of the sinusoidal nature of (|25p . 

Hysteresis in the ring can arise from the self-inductance of the ring as 
follows: The total flux enclosed by the ring will be sum of two terms 

$ = $ H + $ 7 (28) 

where is the flux imposed onto the ring from a pair of Helmholtz coils, 
i.e., by an externally applied H field, and $/ is the flux generated by the 
self-inductance of the ring due to the supercurrent I. For simplicity, we 
shall first consider the case where the width of the ring is less than the 
effective penetration depth, and where we are imposing an external field 
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close to the effective critical field H c , so that the exterior, imposed flux 
can freely flow in or out from the interior of the ring near "criticality" . 

Now the flux generated by the Josephson supercurrent I via the in- 
ductance L of the ring obeys the relationship 

<E>/ = LI (29) 

where L is the self-inductance of the ring. Putting this together with (1251) . 
one obtains the transcendental equation 

/ = /jsin{|^($i f + £/)} (30) 

One can express this equation as two equations in two unknowns, "I> and 
I, as follows: 

7 = /jShl (lH (31) 
$ = $h+LI (32) 

One can rewrite the second equation as the following linear equation for 
/ in terms of <E>: 

7 = I($-$ H ) (33) 

Lj 

which has the form of a straight-line relationship 

I = a$ + b (34) 

where 

a = \ (35) 



(36) 



and 

b = -~ 

are the slope and the intercept of the straight-line relationship 
spectively. 

A graphical method of solution to the transcendental equation ex- 
pressed as the two equations (|31[) and (I34|) is illustrated in Figure 3. 

The intercepts (the black and the white dots) between the two graphs 

/ = Ijsin (|^$) and (37) 

I = a$ + b (38) 

represent the graphical solutions to the transcendental equation. The 
white dot presents an unstable solution (one with an "anti-Lenz" Faraday- 
like law). The black dot solution will eventually settle down to the nearest 
possible metastable, persistent current solution, which is represented by 
the red-dot solution, after the applied magnetic field is reduced to zero, 
and a quantized number of flux lines corresponding to the red-dot solution, 
has been trapped inside the ring. This happens when one reduces the 
external field to values well below the effective critical field H c , where the 
effective penetration depth becomes much smaller than the width of the 
red ring. As one further reduces the applied field until it reaches zero, i.e., 
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7=/ J sin(27iO/(D ) 



I=a®+b 



Figure 3: Graphical method (not to scale) for the solution of the transcendental 
equation ([30|) . written as the two equations (|37]l and (|38|) . The intercept at the 
black dot represents a stable solution, and that at the white dot an unstable 
solution. The red dot represents the nearest possible metastable flux-trapping 
state (i.e., a persistent-current state), which possesses an integer multiple of flux 
quanta $q trapped within the ring, i.e., a nonzero remnant magnetization when 
H — > 0. Since $o is a very small quantity, the sinusoidal oscillations of / are 
very rapid, and the flux-trapping states are highly metastable. 
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Figure 4: SQUID-like model of Figure 2(a) after the insertion of a ferromagnetic 
core (yellow) overcoated by a thick SC coating of lead (blue) into its center (not 
to scale). The flux from the ferromagnetic core contributes to the flux seen by 
the Josephson junction (pink), even though only the vector potential, but not 
the ferromagnetic field from the core, reaches the SC aluminum ring (red). 

as H — > 0, the remnant field as measured by the Hall probe in Figure 1 will 
approximately be the trapped flux corresponding to the red-dot solution, 
divided by the area of the ring. (However, the above analysis is for a 
one-dimensional ring, and needs to be generalized to a three-dimensional 
ring in order to confirm the above conclusions. See Appendix A.) 

The way to extend the above ideas to the configuration in Figure 1 
where there is a ferromagnetic core thickly overcoated by a strong SC 
such as lead or lead-tin solder, is obvious from an inspection of Figure 
4. The above analysis is changed by adding the flux $f c arising from the 
ferromagnetic core to the total flux seen by the Josephson junction, i.e., 



This has the effect of adding a DC bias to the flux in the graphical solu- 
tion of Figure 3 by an amount $fc , thus either increasing or decreasing, 
depending on the sign of $f c , the remnant magnetic field as measured 
by the Hall probe in Figure 1. This analysis also extends to the case of 
a ferromagnetic torus thickly coated with a strong SC (e.g., lead) linked 
with a ring composed of a weak SC (e.g., aluminum or tin), as pictured 
in the first Figure of pQ. 



(39) 
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1 Appendix A: Case of a wide aluminum 
ring. 

An external magnetic field H imposed upon the wide aluminum ring in 
Figure 2 will lead to an inner current /inner flowing along the inner perime- 
ter of the ring within a penetration depth, and an outer current /outer 
flowing along the outer perimeter of the ring within a penetration depth. 
These two currents will join together at the effective Josephson junction 
to satisfy the Josephson sinusoidal relationship 

/ = /inner + /outer = /j sin <f> (40) 

The metastable solutions corresponding to flux quantization in integer 
multiples of the flux quantum, such as the case of the red-dot solution 
in Figure 3, corresponds to solutions with the discrete Aharonov-Bohm 
phase values 

<t> n = 27rn,(n = 0,±l,±2,...) (41) 

i.e., the zero-crossings of the sine function (|40[) . or to the quantized flux 
values 

$„=nA( n =0,±l,±2 r ..) (42) 

The free energy of the system will have extremely deep, local minima 
at these quantized values of the flux, which correspond to the extremely 
metastable solutions for the wide ring at the effective critical values of the 
applied magnetic field (i.e., at "criticality" ) where 

In = /inner | n + /outer \ n = I J sin <f> n = (43) 

In other words, at "criticality", the inner and outer currents will satisfy 
the relationship 

/inner | „ = — /outer |„ (44) 

However, as the externally applied magnetic field H is subsequently re- 
duced to zero, the inner current must remain constant, since the interior 
flux must remain quantized, whilst the outer current will be reduced to 
zero along with the vanishing exterior H field, i.e., 

/inner | „ -> ^ (45) 

as H — > 0, but 

/outerL^O (46) 

as H — > 0. Thus there will result a "remnant" field contained within the 
inner perimeter of the ring, which will approach the value 

-Brcmnant - ► (47) 

as H — > 0, where A is the inner area of ring (excluding the cross-sectional 
area of the torus), which can be measured by the Hall probe (orange) in 
Figure 1. This leads to the possibility of an observation of the "asymmetric 
hysteresis" phenomenon in the linked torus-ring system, as described in 

ra- 
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In summary, the flux measured in the Josephson effect in a SQUID- 
like configuration with a ferromagnetic core is the total flux enclosed by 
this kind of SQUID, including the flux from the ferromagnetic core, as 
illustrated in Figure 4. It is this that leads to the macroscopic, Aharonov- 
Bohm-like effect described here. 
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